Functional minimization method addressed to the vacuum finding 
for an arbitrary driven quantum oscillator 



00 

o 
o 

> 
O 

in 



a,: 

-I— > 

c ■ 

=5 

cr 



> 

m 

On 

in 

o 

oo 
O 



S.V. AnischenkcQ 

Belarus State University, Bobruiskaya 5, Minsk 220050, Belarus 

S.L. CherkasQ 

Institute for Nuclear Problems, Bobruiskaya 11, Minsk 220050, Belarus 

V.L. Kalashniko\@ 

Institut fur Photonik, Technische Universitat Wien, 
Gusshausstrasse 27/387, Vienna A-1040, Austria 
(Received November 5, 2008) 

The old problem exists for a driven (time-dependent) quantum oscillator: to differ the true vacuum 
state from the squeezed one. We suggest finding the true vacuum state by minimization of the 
functional containing the difference of the potential and kinetic energies of oscillator. Analytical 
and numerical examples confirming this offer are considered. 
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I. INTRODUCTION 

A time-dependent (driven) oscillator arises naturally 
in a number of fields of the theoretical physics [J, 
In particular, it has an application in cosmology and 
astrophysics, where the scalar, fermion, gravitational, 
and other quantum fields evolve in an expanding Uni- 
verse 0, 0| . Particle creation by the nonstationary grav- 
itational field is long considered as one of the possible 
sources of the matter origin in the Universe, and, to talk 
about a "particle" one has to understand what is the vac- 
uum. It should also be mentioned that according to the 
modern view the vacuum fluctuations were the seeds for 
the structure formation in Universe 

Nevertheless, the definition of the ground (vacuum) 
state remains to be obscure @, S @, 0|- This forces 
one to use the well-known adiabatic states in concrete 
calculations Q, whereas for systems which admit ana- 
lytical consideration, for instance, quantum field in the 
De Sitter Universe, Bunch-Davis vacuum states Q can 
be built. However it would be desirable to define vac- 
uum state without appealing to the adiabatic series or 
analytical solution (it may be impossible). This issue is 
addressed in our paper. The suggested method allows 
finding numerically the true vacuum state (if it exists). 

Let us remind the problem in more detail. 

Hamiltonian of the time-dependent oscillator has the 
following form: 



H 



1 



f 



u 2 (t)x 2 



(1) 
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Its quantization in the Heisenberg picture consists in a 
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replacement of the coordinate x by the time-dependent 
operator: 



x(t) = ku{t) + k + u*{t), 



(2) 



where the operators a and a obey the commutator re- 
lation 



[a,a+] = I, 

whereas the function u satisfies the equation 



u u — uu = 1. 



(3) 



(4) 



These rules provide the standard commutation relation 
for the momentum and coordinate operators 



x,x 



(5) 



The vacuum state is denned as the state, which is the 
null-space of the annihilation operator: a\0 >= 0. How- 
ever, there remains a problem in the concrete definition 
of u. This function satisfies the oscillator equation of 
motion 



u + u> u = 0, 



(6) 



and one has to define (at some instant) the initial condi- 
tion corresponding to the true vacuum state. It should 
be noted that there exists a family 0, ||| of the func- 
tions u, which satisfy Eq. (U) and are interrelated by the 
Bogolubov's transformation: 

u(i) = cosh r Mo (i) + sinhr e l5 iiQ(t) 
u*{t) =coshru* (t)+smhre- i6 u (t). (7) 

In the terms of an oscillator with constant frequency 
u> = const, this means that it is necessary to differ the 
true vacuum from the squeezed vacuum states. When 



2 



lj — const the vacuum choice can be made by minimiza- 
tion of the mean value of < 0|-ff|0 >, but it is no so for 
time-dependent oscillator. 

However let us remind that when uj — conat the mean 
value of an observable oscillates with time for a squeezed 
vacuum state, whereas it is zero or constant for a true 
vacuum state. One can suppose, that this is the clue to 
issue of the true vacuum state of time-dependent oscilla- 
tor. Namely, an observable value oscillates with time for 
a squeezed vacuum state, but it is monotonic function 
for a true vacuum state. It is convenient to choose the 
difference of the oscillator kinetic and potential energies 
to be this observable. For the vacuum state of oscilla- 
tor with constant frequency, this quantity equals to zero 
according to the virial theorem and we will see that this 
quantity is a monotonic function of time for a vacuum 
state of time-dependent oscillator (if such a state exists) . 



II. TIME-DEPENDENT OSCILLATOR: 
EXAMPLES OF THE VACUUM STATES 
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FIG. 1: The function a 2 for the dependence of u)(t) given by 
((HJ and k = 1, H = 1. Monotonic curve corresponds to the 
vacuum state (r = 0, 8 = ) and oscillating curve corresponds 
to the squeezed state (r = 0.005, 5 = 0). 



We intend to concentrate the different examples, which 
grade with the asymptotic of the adiabatic parameter 
lu/uj 2 : i) uj ~ * const, lu/uj 2 — > 0; ii) uj does not tends to a 
constant but lu/uj 2 — > 0; iii) then lu/lu 2 = const, and, at 
last, iv) uj/lu 2 — > oo. 

Let lu depends on t in the well-known form (2j: 



uj{t) = ky/l + t&nh(Ht). 
The solution of Eqs. (0} and © is 0: 

„ / ik ik . W2k I 



(8) 



2-^1 



i;^ + i; 



1 + e 



2Ht 



, (9) 



where 2-^1(0, b; c; z) is the hypergeometric function [T^]. 
An arbitrary solution from a family of the squeezed vac- 
uum states is given by Eqs. (JT]). 

Mean value of the kinetic and potential energies differ- 
ence is expressed as 

< °\\p 2 - \^x 2 \0 >= hiiu* ~ uj 2 uu*) = &{t). (10) 
Here 

a = i < 0\xp + px\0 >= 7}(uu* +u*u) (11) 

has a sense of additional uncertainty arising in the 
Hciscnbcrg uncertainty relation [ll( : 

<|(p- W ) 2 |x|(i- l0 ) 2 |> >1/4 + ( j 2 , (12) 

where po =< p >, xq =< x > and | > is the arbitrary 
state. For a family of the squeezed vacuum states, in- 
cluding the true vacuum, the inequality (|12p becomes an 
equality. 



Fig. [T] shows the cr(i)-function for different values of 
the parameters r, S. One can see that this function oscil- 
lates at some value of the parameters r, S and the only 
parameter r — results in the monotonic behavior of 
a{t). 

It should be noted that the selection rule have been 
offered Q for a vacuum state as a state having the mini- 
mal uncertainty at each moment of time. As one can see 
from the above example, this rule is not satisfied for the 
vacuum state. Indeed, there exists a region in Fig. [TJ 
where the uncertainty for the slightly squeezed state is 
less than that for the vacuum state. Hence, one has to 
conclude that this selection rule is not valid in the general 
case. 

Our suggestion is to correlate a vacuum state with the 
monotonic time-dependence of the functions a or &. That 
is, for this example, the true vacuum corresponds to the 
function uq given by ©. As the criterium for choosing 
the function with monotonic behavior, one can use the 
minimization of the functional 



Z M )=hm(4^^ 



(13) 



where ro, <5o are some fixed values used for normalization. 
The exact analytic calculation of the functional with the 
function u from |(7J), © and values ro = In 2, Sq = gives 



Z = sinh 2 r. 

225 



(14) 



Thus, the minimization of the functional leads to the 
value r = for the vacuum state. This is because the 
function a has the asymptotic 

a(t) w --sin (2V2kt + S^j sinh(2r) (15) 
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FIG. 2: The function a 2 for -in and -out vacuum states, 
at infinity. 

Instead of the parameters r and S, one can seek the 
initial conditions for Eq. ^ at some t . Really, the rep- 
resention u(t) = e l{p{t) 9{t) leads to (p = -0~ 2 /2 from Eq. 
(]4]). That is, 6*(io) and 0(to) define u(to) and u(io) com- 
pletely because the phase f(to) can be chosen to be zero. 
Then one can solve Eq. ([5]) with some initial conditions 
and find the value of the functional lfT3"|) . Initial condi- 
tions giving the minimum of the functional correspond 
to the vacuum state. 

Moreover, one can write the differential equation di- 
rectly for a(t) [Hj]. 

Straightforward computation shows that a satisfies the 
equation 

a - a - + - + Actio 2 + g %ioio = 0, (16) 

\lO UJ J \ 10 J 

if u(t) obeys 

The relation ((4]) leads to 

— (Actio 2 + cf) Ucrto 3 + if 10 - 2&Co) - 4a 2 = 1, (17) 

for the states belonging to a family of the squeezed vac- 
uum states including the true vacuum. Left hand side of 
Eq. (fT7|) is the integral of motion of (fT6]) . 

One can connect the initial condition for Eq. (|16|) with 
that for Eq. ©: 

a = 66, (18) 
& = 6 2 + 8 2 (<p 2 -lu 2 ). (19) 

Second derivative of a can be expressed through Eq. 
(fT7|) . Thus, the determination of a(to) and &(to) allows 
solving Eq. (fTB|) instead of defining 6(to) and 9(to) and 
solving Eq. ([6]). 

In the above example, the function a has the mono- 
tonic behavior in the vacuum state within all range of 
t. This means that the single global vacuum exists. The 
more complicated case @ is 

= k^2 + t&nh{Ht), (20) 

where there are two different non zero values of w at 
t — > oo and t — > — oo. 
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The dependence u(t) is given by (|20|) and k = 1, H = 1. 

One can see from Fig. [2l that two vacuum solutions ex- 
ist. One of them has the monotonic behavior at t — * +oo 
(out- vacuum state) and the second one has such a behav- 
ior at t — * — oo (in- vacuum state ). In this paper we do 
not discuss an important issue concerning the dependen- 
cies u}(t) providing the unique global vacuum [131 ], but if 
in- or out- vacuums exist, the out- vacuum can be found 
by the minimization of the functional (|13|) , whereas the 
in- vacuum state can be found by setting T — > — oo in 

(USD. 

The considered cases are simple in the sense that u> 
tends to a constant and a notion of particle is asymptot- 
ically defined. For example, the out-vacuum means an 
absence of the particles at t — > oo and, simultaneously, 
the function a has the monotonic behavior. 
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FIG. 3: The function a' 2 for the dependence ui(t) given by 
(|21[) and k = 1, H = 1. Monotonic curve corresponds to the 
vacuum state (r = 0, 8 = ) and oscillating curve corresponds 
to the squeezed state (r = 0.005, 8 = 0). 

Now let us consider the example 

w(t) = VI + Ht, (21) 

where uj(t) does not tend to a constant at infinity, but the 
adiabatic condition 1-^1 — > is still satisfied at t — > oo. 
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FIG. 4: The function a 1 for the dependence co(t) given by 
(|27[) and k = 1, H = 1. Monotonic curve corresponds to the 
vacuum state (r = 0, 5 = ) and oscillating curve corresponds 
to the squeezed state (r = 0.005, 5 = 0). 



Eqs. ©, (|1]) are solvable in the closed form 



i/o 



= H- 1/6 k~ 1/3 V2^Ai 



' 2k 2 ^(Ht + l) \ 
# 2 / 3 (1 - Wf) ) ' 



(22) 



where Ai(z) denotes the Airi function [l(|. Again, the 
calculation of the functional (fl3|) gives Eq. (|T4"|) . It 
should be noted that the asymptotic of er is (see also 
Fig.© 



If (Ak{Ht+lfl 2 _ 
'I C °H ZH + * 



+ V3sin r 4fc ^ + 1)3/2 



3ff 



+ <5 sinh(2r). (23) 



Let us come to the example, where the adiabatic con- 
dition is not fulfilled: 



l + 2ffi 



(24) 



The adiabatic parameter |^-| = 2M. [ s constant. How- 
ever, as it will be shown, the vacuum state exists in this 
case, as well. The solution of Eqs. (O, (|4|) is 

(1 + 2Ht)-' Wk2 - H2 / (2H '> +1 / 2 



V2(k 2 -H 2 y/ 4 
and results in the following asymptotic 



(25) 



2\7fc 2_ 



L===^ffcosh 2 (r) + Hsmh 2 (r) - \J k 2 - H 2 sin^S+^Vk 2 - H 2 hx(2Ht + l)j sinh(2r) 

+iJcos (s+ jj\/k 2 - H 2 ln(2Ft+l)^ sinh(2r)^j . 



(26) 



The comparison with the previous cases demonstrates 
that the constant component appears in the asymptotic 
[Til ]. However, again the functional (fl3"|) has the form 
(fT4")l for k > H . In the opposite case the function a has 
non-oscillating behavior at infinity under arbitrary initial 
conditions. 

Now let us consider the following example: 



k 



1 + HH 2 • 



where the adiabatic parameter | ^ | = becomes 
greater than unity at large t. The expression for the 
function uq has the form 



i/0 



V2<fW 2 



(27) as well as the expression for a is 



; 2 exp = y-^\f H 2 + k 2 arctan(i/t)j , 

(28) 



tcosh(2r)H 2 / H 2 t ( ' 2VH 2 + k 2 arctan(ift) \ 1 /, 2^H 2 + k 2 tan" 1 (iff) \ \ , , , 
' ....iv i _ — sin I H — I I smh(2r) 



2VH 2 + k 2 \ 2VH 2 + k 2 V H 



(29) 



In this example, the function a(t) for an arbitrary state has non-oscillating behavior at infinity (Fig. [3]). It oc- 
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curs because the function u itself ceases to oscillate at 
t — ► oo. In the literature such a phenomenon is in- 
terpreted as the transition from the quantum field to the 
classical one. The behavior of a confirms this interpre- 
tation because the absence of oscillations means the ab- 
sence of interference (i.e., in fact, absence of the main 
constituent of quantum mechanics). In any case we can- 
not talk about an existence of the out- vacuum state here. 
However, one can introduce a concept of the approximate 
vacuum state corresponding to some range of t. One can 
see from Fig. 21 that there exists a range, where the typ- 
ical non-vacuum a oscillates and thus, an approximate 
vacuum state corresponding to the non-oscillating a can 
be defined. 
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FIG. 5: The dependence (|35[1 of the Universe scale factor a(r) 
on the conformal time r. 



III. VACUUMS OF THE SCALAR FIELD 
OSCILLATOR 



In principle a number of the approximate vacuums cor- 
responding to the different time regions can exist. Let us 
take an example, which does not admit an analytical con- 
sideration. 

Lagrangian corresponding to the modes of the scalar 
field in an expanding Universe has the form Q 

Cscai = \:'^a 2 (f)' k (t)'_ k ~a 2 k 2 <t)k(t>-k-a i m 2 (f) k (f)^ k , (30) 



where (f> k is the Fourier-transform of the scalar field 
<M r ) — J2k < fik elkr and a(r) is the scale factor of Uni- 
verse, t is the conformal time dt — a(r)dT. 
The equation of motion can be deduced 



<j>l + (k 2 + a 2 m 2 )cp k + 2-(f>' k = 0. (31) 
a 

Quantization of the scalar field [2] 

^fc = ^-kXl{T) + afeXfc(r) (32) 

leads to the operators of creation and annihilation with 
the commutation rules a£] = 1. The complex func- 
tions Xk{T~) satisfy the relations Q: 



X 'l + (k 2 +m 2 a 2 ) X k + 2-x' k = 0, 
a 

a 2 (r)(Xk x'k -XkXk) = 1 



(33) 



Substitution of Xfc = u k /a results in the time depen- 
dent oscillator: 



u 'k + ( fc2 + fl2m2 H ) u k = 0. 



(34) 



Now we consider some illustrative time-dependence 
a(r) 



a(r) = t (1 + cxp(3 - r)) _1 +4 ( 1 - ' + ' 



1 + cxp(r + 25) 



which is shown in Fig. [51 There exist three ranges, where 
one can search for the vacuum state. Namely, one can 
try to find the in-, out- vacuum states and, besides, the 
approximate vacuum state for the central range (from rs 
—20 to shown in Fig. [5j. The numerical minimization 
of the functional 

Z(a,f3) = f 2 (a T a(r,a,/3)) 2 dr, (36) 

where <t(t) obeys (fToT) . (fTT)) and a = <j{t\), (3 — <t'(ti), 
allows finding the initial conditions corresponding to the 
vacuum. The solutions are shown in Fig. [6l 



IV. VACUUMS OF THE FERMIONIC 
OSCILLATOR 

Let us come to the fermionic oscillator. After decom- 
position of the bispinor ij}(r ) in the complete set of modes 

V'( r ) = J2k i ; k elkr : Lagrangian of the fermion field in the 
expanding Universe (see [H, [lj| [I?], HI] and reference 
therein) takes the form 

L = —iphdrtpk - -^-dri/j^ipk ~ a 3 ip+(ak) ip k 



-a 4 mip k /?Vfc-(37) 

The equation of motion is 
3d' 

iip' k - (ak)ip k + i— ip k - m a(3tp k = 0, (38) 
la 



Fermion field is quantized as 

4>k = bt k>s V-k,s + a k , s u k , s , 
where the bispinor is fl9l ]: 



(39) 



(35) 



Uk,s(v) = 



iXk + ma Xk / 



6 
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FIG. 6: The function a for the vacuum states inspired by the dependence (|35[) and k = 1, m = 1/16. Left column corresponds 
to the scalar field oscillator, whereas right column corresponds to the fermionic one. Upper, middle and bottom rows correspond 
to the solutions representing the vacuums for the out-, in-, and central-time-ranges, respectively. 



and spinors tp s are <y?+ = 



) and v - = ( 1 f 

The bispinor Vk, s is expressed as Vk jS — i7°7 2 (wfc,s) T , 
where the symbol T denotes the transpose vector and 
u = u + 7°. The functions Xk{v) satisfy j 1 81 ] 



,2„2 



Xk + (to + m a — ima 

i *\ 



Xk = 0, (40) 



k 2 XkXl + (arnxl ~ iXfc*) (amXk + ix'k) = h ( 41 ) 

and again there appears time-dependent oscillator (with 
the complex frequency), where the functions Xfc plays a 
role of the above mentioned u k ■ The true vacuum state 
can be defined as that providing a non-oscillating behav- 
ior of the function 



1 



= o (XkXk + Xk Xk) ■ 



(42) 



One can deduce that if Xk{t) obeys (|4H|) then a k satisfy 



/ , 4M"k 2 4M 2 M"\ 
Y1MM' cr fc = 0, (43) 



where M(r) = ma(r). The relation (|4Tj) gives 

(k 2 + M 2 ) (4a k k 2 + Aa k M 2 + a'lf 
-Z^jv'k (^kk 2 + 4a k M 2 + a'i) + 4k 2 a 2 + a' 2 = 1. 



(44) 



The vacuum solutions obtained by minimization of the 
functional (|36f> in the three different ranges are shown in 
Fig. [5] (right column). 
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V. VACUUMS OF THE TWO COUPLED 
OSCILLATORS WITH CONSTRAINT 

Now we address ourself to a little more complicated 
system: namely, the two time-dependent coupled oscilla- 
tors with constraint. This system appears in the theory of 
anisotropy of the Cosmic Microwave Background [22L I23I] . 
One can expect, that some difficulties will arise with the 
vacuum definition for this system, because the quantiza- 
tion of constrained systems can reveals some nontrivial 
features. However, we will see that there are no patholo- 
gies in this particular case. 

Both scalar field and gravitation can be assumed to be 
specified by the action [2, Q 



S = - 



1 



]( ._ r , ; &\V~gR+ f dW^K^) 2 -^)], 

45) 



Representation of the metric tensor in the form 0, H 

ds 2 = (1 + 2$(r, t))dt 2 - a 2 {t)(l - 2$(r, t))dr 2 , (46) 

and considering the scalar field as that possessing a 
spatially uniform component with a small perturbation 
around it: 



{r,t) = <t>{t)+e{r,t) 



(47) 



allows obtaining the system of equations |3[ of zero order 
in 9 and $, 



- a z a + </>V + 2a 6 V{4>) = 0, 



a = ~-a<p - — + 3aV((f>), 
2 2a 

+ 3^+^ = 0, 
a d(f> 



(48) 
(49) 

(50) 



where we use system of units AttG/3 = 1. The first 
order equations for the Fourier-transformed perturba- 
tions of the scalar field 9(r, t) = 2~2 k ®k(t)e lkr and metric 
®(r,t) = J2k®k(t)e ikr have the following form 0,112]: 



\®kk 2 + ^-0 k a 2 + $ fc aa + 9 k ^a 2 + 2<S> k a 2 V {4>) = 0, 
3 dq> 

(51) 

~*fc-^+M = 0, (52) 
3 3a 

- 3 lT fe - ** - 4$ fc - + 30 fc - 6V(0)$ fc = 0, (53) 
acp a 

a k 2 d 2 V dV 
k + 3-8 k + —9 k + — 9 k + 2— $ fc - 40$ fc = 0. 
a ar a<p z dip 



(54) 



Eqs. (|53|. ([54]) are the equations of motion. They can 
also be obtained from Lagrangian 



1 d "V d"V 1 

L = ~ oTHa fl * fl -*° 3 + 2— ^$- fc a 3 - 10V$ fe $_ fc a 3 - -ak 2 k 9_ k 
A — ' I &(b A ad) 1 

k 

-ia/s 2 $ fc $_ fc + ^-a 3 9 k 9- k - \a 3 $ k <S>- k - 4$ fc $_ fc a 2 a - 4a 3 $ fc 0_ fc 0. (55) 
O 2 2 



Eqs. (|51| . (f52|) are the constraints. However, Eq. ([51]) 
can be derived from Eqs. (f52")) . (f53")) . ([54")) and, thus, it 
is not independent. That is there are two time depen- 
dent oscillators with one constraint of the first kind [2J] . 
Using this constraint one can exclude the scalar field per- 
turbation from (|53p and obtain 



*4t 



a d 
a dt 



k 2 


2a 


2d 2 


^ + 


a 


~a^ 














(- 




U 2 


a J 





where the uniform scalar field <p has been excluded by 
using l|48[). (j49|) as well. Quantization consists in 



-k u l(t) +&kU k {t), 



(57) 



and the vacuum can be found by minimization of the 
quantity Z(a k ,(3 k ) (|36p . that allows finding a k and (3 kl 
which correspond to the vacuum state. In this we have 
solved Eq. (|56[) directly and have written the initial con- 
ditions at t as Uk(to) = a*, «fc(*o) = Pk~ i^T^T' wnere 
■X"(*o) — ^3 ~ "S-L . • These initial conditions are con 



a J a 2 lt=t 

sistent with the relation 



iX(t), 



(58) 



WfcU fc - u k u k 

which is analog of (|33)) and corresponds to the general 
case of quantization of oscillator with the time-dependent 
mass and frequency [251 ]. 

On the other hand, one can express through 9k, 
and obtain the equation for 9 k analogously to (|5"6j) . This 
equation turns out to be more complicated and we do 
not write it here. The question arises: would be the 
vacuum state the same, if the scalar field perturbation is 
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FIG. 7: The function a 2 corresponding to the vacuum for 
a(t) = sinh(7t), where 7 = and k = 1. top panel shows 
result of the metric perturbation quantization, bottom panel 
shows the result of the scalar field perturbation quantization. 

quantized as: 

e k = &t k U* k {t) + & k U k {t)l (59) 
From Eq. (|52p . we have 

U k = — r ( u k + -u k ) , (60) 
30 V a J 



where Eqs. I|48[). (|49|) reduce to the form = 
/1 _ is 

V 3 a 2 3 a' 

Let us consider the particular case of a(t) = sinh(Tt). 
The numerical minimization of the functional allows find- 
ing the vacuum solution for which the function u\ k — 
\{u k u* k + u k u k ) has monotonic behavior and, thus, cor- 
responds to the vacuum state. Now if one expresses 
(J 2k = \(UkUl + U k U k ) through the functions Uk given 
by (jfjfj)) . it is seen from FigJTl that (<J2k) 2 has also mono- 
tonic behavior. Thus, it is the vacuum state for the k - 
oscillator too. Moreover, one can choose any convenient 
variable from a combination of $ k and 9 k as it usually 
done HiHl lift 

Procedure " NMinimize" of the Wolfram software 
" Mathcmatica" is used in all the numerical calculations. 



VI. CONCLUSION 

We have considered the method to find the vacuum 
state of a driven quantum oscillator numerically by the 
means of minimization of the functional containing the 
square of derivative of the additional uncertainty a aris- 
ing in the Heisenberg uncertainty relation. For a time- 
dependent oscillator, the derivative of a coincides with 
the difference of kinetic and potential energies. We show 
that this method can also be applied to both fermionic 
oscillator and pair of the coupled constrained oscillators. 
The last example is widely used in the theory of the mi- 
crowave background anisotropy. We have verified that 
there is no problem with a selection of the vacuum state 
for the last system in spite of some discussions appearing 
in the literature [13, HI] • 
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